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ABSTRACT 


V* 


A  mathematical  model  of  the  range-gate  pulloff  electronic  counter¬ 
measure  is  developed  based  on  a  statistical  model  for  the  radar  return  from  a 
slowly  fluctuating  point  target  in  the  presence  of  white  noise.  It  is  shown 
that  the  same  model  is  also  appropriate  for  describing  velocity-gate  pulloff. 
An  optimization  problem  is  formulated  which  determines  the  jammer  delay 
inducing  the  maximum  bias  in  the  range  estimation  processor  of  the  victim 
radar.  These  results  are  then  applied  to  the  specific  case  in  which  the 
transmitted  signal  is  a  gaussian  pulse.  The  optimal  delay  and  bias  are 
calculated  as  functions  of  the  signal-to- jammer  power  ratio  and  the  pulse 
width  of  the  transmitted  signal. 


RESUME 


Un  module  mathgmatique  de  contre-mesure  glectronique  par  rgtabllssement 
des  portes  distance  est  Slaborg  A  partir  d'un  mod&le  statistique  applicable 
aux  echos  radars  provenant  d'une  cible  ponctuelle  qui  fluctue  lentement  en 
presence  de  bruit  blanc.  II  est  dgmontrg  que  le  m§me  mod&le  peut  tlgalement 
decrlre  le  retablissement  des  portes  vitesse.  Un  probl&me  d' optimisation  est 
formula  pour  determiner  le  retard  du  brouilleur  introduisant  un  biais  maximal 
dans  le  calcul  de  la  distance  par  le  processeur  du  radar  victime.  Les 
rgsultats  sont  ensuite  appliques  au  cas  particulier  oO  le  signal  transmis  est 
une  Impulsion  gaussienne.  Le  delai  maximum  et  le  biais  sont  calculus  en 
fonction  du  rapport  des  puissances  signal/brouilleur  et  de  la  dur£e 
d' impulsion  du  signal  transmis. 


(ill) 


-.v v.v.v v./. 


TABLE  OF  CONTENTS 


ABSTRACT/RESUME . iii 


TABLE  OF  CONTENTS. 


LIST  OF  FIGURES . vii 


1.0  INTRODUCTION. 


2.0  BASIC  RADAR  MODEL 


3.0  VELOCITY-GATE  PULLOFF 


4.0  DETERMINING  OPTIMAL  BIAS. 


5.0  CONCLUSIONS . 21 


6.0  REFERENCES . 22 


APPENDIX  A:  COMPUTER  PROGRAMS  .  A-l 


Tf  H  »  (T»  V 


LIST  OF  FIGURES 


FIGURE  1 


FIGURE  2 


FIGURE  3 


FIGURE  4 


FIGURE  5 


GRAPH  OF  FUNCTION  g1 .  12 


GRAPH  OF  FUNCTION  g3. 


COMPARISON  BETWEEN  Tr,  T  AND  T*. 


GRAPH  OF  g4 


GRAPH  OF  g3  FOR  SEVERAL  VALUES  OF  Tr/Ts  SUPERIMPOSED 

ON  GRAPH  OF  g! .  19 


FIGURE  6: 


GRAPH  OF  SOLUTIONS  Topt  AND  topt, 


(vii) 


*  ** 


l 


1.0  INTRODUCTION 


DREO  has  a  continuing  interest  in  simulating  hypothetical  engagements 
between  fighter  aircraft  and  one  or  more  surface-to-air  missiles.  One  of  the 
aims  in  developing  such  simulations  might  be  to  discover  whether  overall 
performance  of  the  aircraft  can  be  enhanced  by  providing  it  with  an  electronic 
countermeasures  (ECM)  suite.  In  particular,  two  of  the  countermeasures  which  are 
frequently  encountered  in  the  open  literature  [1.]  are  range-gate  and 
velocity-gate  pull- off.  In  this  technical  note  we  will  propose  a  simple  model 
which  can  be  used  to  describe  either  of  these  phenomena.  We  will  then  use  this 
model  to  derive  some  general  "rules  of  thumb”  regarding  the  possible  performance 
and  response  time  of  a  hypothetical  ECM  system  which  uses  these  countermeasures. 

2.0  BASIC  RADAR  MODEL 


We  shall  first  outline  a  statistical  model  for  a  slowly  fluctuating  point 
target  in  the  case  of  white  bandpass  noise.  Since  this  model  is  described  in 
Chapter  9  of  [2]  in  some  detail,  we  shall  dispense  with  derivations  and  formal 
proofs,  and  restrict  our  presentation  to  a  summary  of  the  principal  assumptions 
and  results.  In  particular,  we  will  make  the  following  basic  assumptions: 
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(i)  the  aircraft  can  be  modelled  as  a  number  of  reflecting  surfaces, 
such  that  the  return  from  each  of  these  can  be  described  as  the 
product  of  an  independent  complex  gaussian  random  variable  and  a 
time-varying  (complex-valued)  sinusoid.  Moreover,  we  assume  that 
there  are  sufficiently  many  reflecting  surfaces  that  the  central 
limit  theorem  can  be  applied  to  the  sum  of  reflected  signals; 

(ii)  the  reflection  process  is  linear  and  frequency-independent; 

(iii)  the  change  in  time  scale  of  the  complex  envelope  due  to  non-zero 
target  velocity  can  be  ignored  (i.e.  only  the  Doppler  shift  needs  to 
be  considered); 

(iv)  all  other  stochastic  variations  in  the  radar  returns  can  be  modelled 
as  additive  white  bandpass  noise. 

Subject  to  these  assumptions,  we  can  derive  the  following  mathematical 
description  of  the  received  waveform  In  the  presence  of  a  moving  target: 

s(t)  =  Re  [r(t)  exp  (jwct)]  (2.1) 

where  r(t)  is  defined  by 
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r  (t )  = 


bs  /Es 


f(t-  us)  exp  (jws  t)  +  n(t ) 


(2.2) 
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and  the  following  notation  and  conventions  ai j  adopted: 


(1)  f(t)  is  the  complex  envelope  of  the  transmitted  signal 

(2)  wc  is  the  carrier  frequency 

(3)  us  is  a  delay  proportional  to  target  range 

(A)  ws  is  Doppler  shift  due  to  relative  target  velocity 

(5)  Es  is  the  transmitted  energy 

(6)  bs  is  a  zero-mean  complex  gaussian  random  variable 

(7)  E(bsbg*)  =  2  v|  (where  the  value  of  Vs  depends  on 

antenna  gains,  path  losses,  target  radar  cross-section,  etc., 
denotes  complex  conjugate,  and  E(  )  indicates  expected  value). 

(8)  n(t)  is  an  independent  zero-mean  white  complex  gaussian  random 
process 

and 

(9)  E(n( t)  n* ( t) )  -  N0  5  (t-s) 
where  we  use  the  notation 

n*(t)  ■  (n( t) )* . 

Clearly  we  can  work  with  either  (2.1)  or  (2.2);  for  convenience,  we  will 
usually  use  (2.2)  in  the  remainder  of  this  paper. 

We  shall  restrict  our  attention  to  the  case  of  transmission  of  a  single 
pulse.  For  this  case,  ug  and  ws  may  be  considered  as  constants;  for 
algebraic  simplicity,  we  can  in  fact  assume  without  loss  of  generality  that 

us  =  ws  =  0  (2.3) 

Hence  we  can  rewrite  (2.2)  as 

r(t)  =  bs  .'Eg  f(t)  +  n(t)  (2.4) 

According  to  [1],  Range-Gate  Pulloff  is  a  self-screening  ECM  technique 
for  use  against  pulsed,  noncoherent,  automatic  range-tracking  radars.  The 
victim  radar's  signal  is  received,  amplified,  and  then  retransmitted  with  a 
minimum  delay,  in  an  attempt  to  provide  a  strong  beacon  signal  which 
"captures"  the  victim's  range-gate.  The  time  delay  in  the  repeated  signal  is 


then  successively  increased  on  a  pulse-by-pulse  basis,  creating  a  series  of 
false  targets.  When  the  victim  range-gate  has  been  moved  sufficiently  far 
away  from  the  true  target  position,  the  ECM  repeater  is  turned  off,  and  the 
victim  radar  will  break  its  range  track. 

In  the  context  of  our  radar  model  (2.4),  we  can  describe  the  deception 
pulse  by  adding  an  extra  term  consisting  of  a  delayed  and  amplified  version  of 
the  complex  envelope  f.  Since  the  deception  pulse  will  in  general  be  out  of 
phase  with  the  transmitted  signal,  we  will  also  include  a  random  phase 
factor.  Indeed,  without  increasing  the  computational  complexity,  we  can 
generalize  the  random  phase  component  to  be  a  complex  gaussian  random 
variable.  Hence  we  obtain. 


r(t)  =  bs  .''Is  f(t)  +  bd  f(t-T)  +  n(t)  (2.5) 

where  T  is  the  delay  in  transmitting  the  deception  pulse,  and  bd  is  an 
independent,  zero-mean  complex  gaussian  random  variable  such  that 

E(bdbd*)  =  2V|  (2.6) 

and 

E(bdbd)  =  E(bsbs)  =  E(bsbd*)  =  E(bdbs*)  -  0  (2.7) 

We  can  model  the  action  of  the  radar  receiver  in  estimating  the  true 
value  of  the  return  signal  as  a  time-invariant  linear  filter  with  impulse 
response  h(t).  For  the  moment,  we  will  only  assume  that  h  is  an  L^(R) 
function  (i.e.  is  square- integratable) .  We  will  be  more  specific  as  to  other 
properties  of  h  later.  Denoting  the  receiver  output  as  r(t),  we  obtain  from 
(2.5) 

r(t)  =  bs  vCsf(t)  +  bd  /Td  f(t-T)  +  n(t)  (2.8) 

where 

f(t)  =  (f*h)(t)  =  r  f ( u)  h(t-u)  du 

n(t)  =  (n*h)  ( t) 

Whether  ***"  denotes  complex  conjugate  or  convolution  will  usually  be  obvious. 
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Now,  we  define  th<_  expected  power  envelope  of  the  filtered  return  signal 

to  be 


P(t;T)  =  E(r(t)  f*(t))  =  (bs  /  Es  f(t)  +  bd  /Td  f(t-T)  +  n(t)) 

x  (bg*^  f(t)  +  bd  >¥d  f(t-T)  +  n(t))* 

=  2 V|  Es  f(t)f*(t)  +  2v£  Ed  f(t-T)f*(t-T) 

+  E(n( t)n*(t) )  (2.9) 


This  last  follows  from  (2.6),  (2.7),  and  assumption  (7)  above.  We  now  wish  to 
evaluate  the  last  term  in  (2.9): 


E(n(t)n*(t))  =  E  [;'  n(u)  h(  t-u)  du][/  n(s)  h(  t-s)  ds]* 

R  R 

=  -r  f  h(t-u)  h*(t-s)  E(n(u)n*(s))  duds 


-  J  h(t-s)  f  N0  6(u-s)  h*(t-u)  du  ds 

R  R  (assumption  (9)) 


=  N0  /  h(t-s)  h*(t-s)  ds 

R 

=  N0  f  h(s)  h*(s)  ds  (translation  invariance  of 
R  Lebesque  integral) 

-  N0  ||  h  | | \  (2.10) 


Substituting  (2.10)  into  (2.9),  we  obtain 

P(t;T)  =  2V|  Es  f(t)f*(t)  +  2v£  Ed  f( t-T)f*( t-T)  +  N  ||  h|(' 

(2.11) 

As  is  pointed  out  in  [3],  we  can  model  the  estimation  of  range  by  the 
tracking  radar  as  a  simple  optimization  problem.  In  effect,  we  claim  that 
those  values  of  t  at  which  P(t;T)  attains  local  maximal  will  correspond  to 
estimates  of  target  ranges.  This  models  leading  edge  tracking  In  the  sense 
that  it  analyzes  centroid  biasing  caused  by  the  addition  of  the  deception 
pulse.  We  note  in  passing  that  this  algorithm  would  not  accurately  model 
tracking  gates  with  sophisticated  statistical  signal  processing  or  threshold 
as  well  as  differencing  logic. 
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We  have  not  thus  far  committed  ourselves  to  a  choice  of  the  receiver 
Impulse  response  h.  If  we  wished  to  extend  our  analysis  to  include  a  more 
complex  estimation  model,  then  a  matched  filter  or  a  filter  based  on  estimates 
of  the  statistical  properties  of  the  interfering  signals  might  be 
appropriate.  However,  given  that  we  are  using  the  simple  estimation  model 
described  above,  we  will  assume  that  the  impulse  response  of  the  radar 
receiver  is  (nearly)  a  delta  function.  In  effect,  this  corresponds  to  a 
receiver  which  estimates  the  target  range  based  only  on  the  return  signal 
without  doing  any  signal  processing. 

Taking  h  as  a  delta  function  appears  to  conflict  with  our  earlier  assumption 
that  h  L2(R).  However,  we  will  now  show  that  there  is  a  sequence  of  L2(R) 
functions,  denoted  d^(t),  which  "look  like"  delta  functions  for  sufficiently 
small  values  of  h. 


Let  d^(t)  be  defined  by 


dh(t)  =  J 


*  & 


(2.12) 


where  a(t)  is  any  non-negative  bounded  continuous  function  with  compact 
support  contained  in  the  interval  [-1,1].  Assume  that  the  total  mass  of  a  is 
1.  Then  following  the  same  style  of  argument  as  is  used  in  (4),  Chapter  6,  we 
can  show  that  for  any  function  frL2(R),  (d^  *f)  converges  to  f  both 
pointwise  and  uniform'’  in  L2(R)  as  h  >0.  For  consider: 


(dh  *f)(t)  -  f(t)P  = 


!  /  -  a  (t~) f ( t-u)  du  -  r  f(t)  a(u)  du  I  ?- 

R  h  h  R 

=  |  /  a(s)  f(t-hs)  ds  -  /  f(t)  a(s)  ds I ' 

R  R 

(using  the  change  of  variable  s  =  u/h) 


[  ■  f(t-hs)  -  f(t)  a(s)  ds]2 

R 

[  f  |  f  ( t-hs)  -  f(t)]2  a(s)  ds][|{  a(s)  ds] 
R  (by  the  Schwarz  inequality) 


(2.13) 


j  a|  1  [.'  jf(t-hs)  -  f(t)|2  ds] 


!  a  ; |  w2  ( f , h) 


(2.14) 


2  ^  •»  v  v 


where  w(f,h)  is  the  oscillation  of  f c  L2  (R) ,  defined  by 

w(f,h)  =  sup  (  /  |f(t-s)  -  f  C  s)  ]  2  ds)-^2 

|s |  1  h  R 

and;  tal  !*,  denotes  the  maximum  value  of  a.  We  note  from  [4]  that  w(f,h)+0  as 
h>0.  This  proves  that  (d^  *f)  converges  pointwise.  Now,  using  (2.13)  and 
integrating  both  sides  we  obtain 

/  | (dft  *f)(t)  -  f(t)  |2  dt  =  f  f  |f(t-hs)  -  f(t)j  2  a(s)  dtds 

R  R  R 

1  w2(f,h)  I  a(s)  ds 

R 

(using  definition  of  oscillation) 

I  w2(f,h)  (2.15) 

which  proves  convergence  in  L2  (R) . 

We  are  now  in  a  position  to  state  the  problem  we  wish  to  consider:  given 
the  model  as  described,  choose  the  value  of  the  deception  delay  T  which 
achieves  the  maximum  range  deception  while  maintaining  (statistically) 
non-resolvable  signal  and  deception  pulses.  Intuitively  this  problem  will 
admit  one  or  more  solutions.  Since  we  can  assume  that  the  signal  envelope 
f(t)  has  been  designed  to  have  only  one  maximum,  very  small  values  of  T  will 
still  generate  only  one  maximum.  On  the  other  hand,  very  large  values  of  T 
will  certainly  result  in  a  power  envelope  which  resolves  the  return  signal  and 
deception  pulse  as  two  distinct  targets.  Clearly,  there  must  be  some  "happy 
medium" . 

According  to  our  model,  the  estimated  target  range(s)  will  correspond  to 
the  maximum  or  maximal  of  the  function  P(t;T)  defined  in  (2.11).  Since  the 
last  term  in  (2.11)  is  constant  (depending  only  on  the  choice  of  the  impulse 
response  of  the  receiver  filter  and  the  spectral  height  of  the  noise),  we  can 
effectively  disregard  it.  Moreover,  we  showed  above  that  we  can  choose  the 
impulse  response  of  our  receiver  filter  so  as  to  make  the  function  f(t)  as 
close  as  we  like  (pointwise)  to  the  function  f(t),  by  choosing  the  parameter  h 
sufficiently  small.  This  implies  that  the  maximum  of  the  function  P(t;T) 
defined  by 

P(t;T)  =  2V2  Es  f(t)f*(t)  +  2V^  Ed  f ( t-T)f*( t-T)  (2.16) 

will  also  maximize  P(t;T).  It  Is  this  function  we  will  work  with  in  section 

4.0. 
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3 . 0  VELOCITY-GATE  PULLOFF 

As  indicated  in  the  introduction,  we  wish  to  extend  our  analysis  to 
include  the  velocity-gate  pulloff  ECM  as  well.  As  described  in  (1),  this  is 
another  self-screening  ECM  technique.  The  signal  from  the  victim  radar  is 
received,  amplified  coherently,  and  retransmitted  to  provide  a  strong  beacon 
which  captures  the  velocity-gate  of  the  radar.  The  Doppler  frequency  of  the 
deception  signal  is  then  moved  away  from  the  true  target  Doppler  frequency  at 
a  rate  that  does  not  exceed  the  victim  radar's  tracking  capability.  At  some 
point  the  ECM  repeater  is  turned  off,  causing  the  radar  to  break  track.  It 
should  be  noted  that  if  velocity-gate  pulloff  is  used  in  conjunction  with 
range-gate  pulloff,  the  rate  at  whic  the  false  range  is  changed  must  equal 
the  false  velf  :ity.  Although  it  may  be  a  subject  of  future  research,  the 
analysis  presented  in  this  paper  cannot  be  applied  to  simultaneous  use  of  the 
two  countermeasures. 

Referring  to  the  basic  equations  (2.2)  and  (2.4)  of  the  last  section,  we 
see  that  a  possible  model  for  velocity-gate  deception  is  given  by 

r(t)  =  bs  ,  Eg  f(t)  +  b^  ,'E^  f(t)  exp  ( jw^  t)  +  n(t)  (3.1) 

under  the  same  assumptions  as  were  posed  in  the  last  section. 

Now,  we  will  denote  the  Fourier  transform  of  a  function  f(t)  by  F(jw)  or 
by  F (f).  Recalling  that 

f'(f(t)  exp(jw^t))  =  F(jw  -  jw^) 
we  obtain  (formally)  from  (3.1) 

R  ( jw)  =  bs  .  Es  F  ( jw)  +  b^  .  F.^  F(jw  -  jwd)  +  N(  jw)  (3.2) 

Now,  from  (3.1)  and  Plancherel's  Theorem,  the  energy  in  the  return  signal  is 

‘ E( r( t) r*( t) )  dt  =  E(R( jw)R*( jw) )  dw  (3.3) 

k 

Of  course,  the  left  hand  side  of  (3.3)  will  not  converge;  however,  using 
an  argument  similar  to  that  of  the  previous  section,  we  can  overcome  this 
technical  difficulty.  Hence,  if  we  define 


S(w)  =  E(R( jw)R*( jw) ) 


=  2v|  Eg  F(jw)F*(jw)  +  2V§  Ed  F( j(w-wd) )F*( j(w-wd) )  +  N0 


(3.4) 


then  we  can  interpret  S(w)  as  the  expected  spectral  power  density  of  the 
return  signal.  Following  the  same  line  as  was  taken  in  the  last  section,  we 
now  claim  that  those  values  of  w  at  which  S(w)  attains  local  maximal  will 
correspond  to  estimates  of  target  velocity.  Since  the  constant  term  in  (3.4) 
due  to  the  spectral  height  of  the  additive  noise  will  not  affect  the  position 
of  any  maximum  values,  we  can  effectively  ignore  it.  Consequently,  the  form 
of  (3.4)  is  identical  to  that  of  (2.16). 

In  the  next  section  we  will  show  how  our  model  can  be  used  to  determine 
the  optimal  bias  for  the  range-gate  problem,  i.e.  equation  (2.16).  However, 
as  we  have  shown,  the  results  may  equally  be  applied  to  the  velocity-gate 
problem,  although  we  shall  not  explicitly  do  so  in  this  paper. 


4.0  DETERMINING  OPTIMAL  BIAS 


In  order  to  illustrate  the  utility  of  some  of  these  ideas,  let  us  assume 
that  the  transmitted  signal  is  a  gaussian  pulse,  i.e. 


- ,  v  ,  l  N 1/4  ,-6“., 

t  ( c)  =  (~pr)  exp  (yp-) 

s  s 


Thus  from  (2.16)  we  obtain 


(4.1) 


P(t;T)  =  2V2  Es  (^7  )1/2  exp(-t2)  +  2vg  Ed  (“r)1/2  exp(  -  (t-T)2) 

s  j  T~  '  s 

s 

=  Es  exp(-t2)  +  Ed  exp(-(t-T)2) 


(4.2) 


Now,  for  each  value  of  T,  the  value(s)  of  t  which  maximize  P(t;T)  must  satisfy 


_  P(t;T)  =  _r2t_ Es  exp(-tz)  -  2(t-T)  Ed  exp(-(t-T)2) 

4-  t2  t2  t2 


(4.3) 


Simplifying  (4.3)  we  obtain 


t  =  exp(t2  -  (t-T)2)  =  expf  1  (2tT  -  T2)  } 


(4.4) 


-V  -\ 


-VVV"-  '*■  '■  *"*  *>  *■«  *'*«  v\*'.  *’  -  «  ",  **-  *  -  v 
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where 


(Es/Ed)  = 


We  note  that  since  the  right-hand  side  of  (4.4)  is  always  strictly 
positive,  and  since  we  can  assume  without  loss  of  generality  that  t  0  ,  we 
must  have 


T-t  0  or  t  <  T 


(4.5) 


Let  us  designate  the  value  of  t  at  which  a  local  maximum  occurs  by  t*. 
Then  in  order  to  guarantee  that  t*  satisfying  (4.4)  jj;  a  local  maximum  (rather 
than  a  local  minimum),  we  impose  the  further  necessary  condition 


9 - P  (t;T)  I  •  0 

~'t-  t  =  t* 


•-  exp  (-t  [-E  +  2E^  t*~  -  E,  €>xp  .  t* ~  -  ( tft-T)  “ , 

T-  T-  s  T  ‘  'I'.:  T.  1 

s  s  s  s  s 


(4.6) 


+  2E  j  (  f-'-T)  “  exp  r  t  "  -  ( t  -T )  “  ,  ]  •  0 


Substituting  (4.4)  into  (4.6)  and  simplying,  we  obtain 


2Eb  exp(-t*2)  [-1  +  2t*2  -t*  +  2(T-t*) t* ] 

Tg  t|  T-t*  T2 

2ES  exp(-t*2)  t  -T  +  2Tt* ]  <  0 

T|  “Tf  T-tA  tT 


(4.7) 


Since  the  leading  factor  in  (4.7)  is  already  strictly  positive,  it  is 
clear  that  the  sign  of  the  left-hand  side  in  (4.7)  depends  only  on  the 
expression  in  square  brackets.  Moreover,  since  T  is  a  delay,  we  can  also 
assume  that  T  0,  resulting  in  the  following  condition  which  is  equivalent  to 
(4.7): 


-1  +  2t*  <  0 

T-t*  TT 


(4.8) 


IT  *  *  *  TJ'f'l  *  *  V*  U  f  'P'»  TVTT*  V*  V  'f  »J  r*J  fw  TV 
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2t*2  -  2t*T  +  T2  >  0 


(4.9) 


where  we  have  used  (4.5)  in  (4.8)  and  simplified.  We  note  that  for 
(T/Ts)  <  /"T,  then  (4.9)  will  be  true  for  all  t*c[0,T].  Hence,  for  all 
choices  of  T  such  that  (T/Ts) <  /2,  we  are  quaranteed  there  will  be  a  single 
local  maximum,  hence  t*  will  designate  a  global  maximum.  We  can  interpret 
this  as  meaning  that  only  one  target  will  be  resolved  by  the  victim  radar. 
However,  it  is  unclear  how  much  larger  T  can  be  chosen  such  that  the  radar 
will  continue  to  resolve  only  one  target.  This  obviously  will  depend  on  the 
actual  value(s)  of  t*  which  solve  (4.4),  which  will  in  turn  depend  on  the 
"power  ratio”  (Eg/E^) ,  and  on  the  "pulse  width"  Tg. 

Now,  in  order  to  resolve  two  peaks,  there  must  be  some  value  of  T  which 
yields  a  "flat  spot",  i.e.,  a  pair  (t*,T)  such  that 

;.P(t*,T)  =  j2P(t*,T)  =  0 
~  Tt2 


From  (4.4)  and  (4.9),  this  gives  the  pair  of  equations, 

2t*2  -  2t*T  +  T2  =  0  (4.10) 

and 

t*  =  exp  f  1  (2t*T  -  T2),  (4.11) 

1  T-t* 


To  simplify  the  solution  of  this  pair  of  equations,  we  introduce  a 
"trick"  based  on  (4.5).  Since  we  know  that  t*<T,  there  is  a  number  F 
(possibly  depending  on  T)  such  that 

0  -  -  1  and  ”T  *  t*  (4.12) 


Substituting  (4.12)  into  (4.10)  we  obtain 


2  2T2  -  2.T2  +  T2  =0 

(1-  )  =  T2  (4.13) 

21  * 


V.K 

s. 

s.'s; 

•* 


Substituting  (4.12)  into  (4.11)  yields 


=  exp  {  (2T2/Tg)  (  5'  -  1)} 


(4.14) 


Now  combining  (4.13)  and  (4.14)  gives 


-■»  =  1  ~  ■'  exp  (  F -  h  )  -  gi(" ) 

F(l-F) 


(4.15) 


Hence  for  a  given  value  of  a  ,  we  can  solve  for  F  =  gj^  (a)  in 
(4.15),  and  then  use  (4.13)  to  find  the  value  of  T,  say  Tr(a),  at  which 
targets  will  be  resolved.  Figure  1  shows  a  graph  of  g^;  we  note  that 
gj-*-  (c*  )  is  defined  and  one-to-one  for  all  a  >  0. 

However,  what  we  are  seeking  is  not  necessarily  Tr(ct  )  ,  but  rather  the 
value  of  T  (less  than  or  equal  to  Tr  (a))  which  gives  rise  to  a  single 
target  with  the  largest  "bias"  value  t*.  In  order  to  find  a  necessary 
condition  which  defines  these  values  of  T  and  t*,  let  us  differentiate  (4.4) 
implicitly  with  respect  to  T: 


t  *  -  qt*  (1-t')  =  1  (2t'T  +  2t*  -  2T)  exp  {  1  (2t*T  -T2)}  (4.16) 

a  T-t*  (T-t*!2  t|  t| 

where  t  =  _d  t*(T).  Substituting  (4.4)  into  (4.16)  yields 
dT 


a  rT  -  at* 

(T-t*)2 
t  T  -  2Tt* 


=  _2_r  t'T  +  t* 

T2  { 

J-s 

(T-t*)}  =  t*{l 


-  T 


H 


a  t* 
T-t* 


} 


2  (T-t*) 

t2 

l8 


(4.17) 


Hence,  we  deduce  that 


_d_  t*(T)  -  t* (1  -  (2/T2)(T-t*)2  I  (4.18) 

dT  [T  -  2  t  t*(T-t* )  ] 

jZ 

s 

To  find  a  (local)  maximum,  we  let  d_  t*(T)  =  0  in  (4.18),  obtaining 

dT 

t*  =  0  or  1  -  (2/T2)(T-t*)2  =  0 


13  - 


Since  we  have  assumed  that  t*  >  0,  we  are  left  with 


T  -  t*  =  Ts//2 


Now,  combining  (4.19)  with  the  basic  equation  (4.4)  yields 


(4.19) 


a  v  2  t* 
To 


{ (t*/Ts)2  -  1  } 


a  =  1  (Ts)  exp  {  (t*/Ts)2  -  1}  =  g?(t») 

*^2  — UF  2  To 


(4.20) 


For  any  value  of  a  >  1,  there  are  two  numbers  t*  such  that 
t*  =  gj}-  (a),  as  shown  in  Figure  2.  For  each  of  these  (4.19)  can  be  used 
to  calculate  corresponding  values  of  T,  say  T*(a)  and  T*(«),  such  that  T*<T*. 

However,  as  shown  in  Figure  3,  T*(,t)  is  actually  greater  than  Tr(rt), 
i.e.,  if  T  =  T*(ot  )  then  two  targets  will  be  resolved.  Hence  we  only  use  that 
portion  of  g2  to  the  left  of  the  dotted  line  in  Figure  2  to  find  t*. 

If u  <  1,  then  by  inspection  there  is  no  solution  for  (4.20),  and  hence  no 
value  of  T  such  that  d_  t*(T)  =  0.  Referring  to  (4.18),  if  we  assume  that 
T  <  l_Tg,  we  have  dT 


t*{l  -  (  2  )  (T-t*)2}  >  t*fl  -  (  2  )  T2 


(4.21) 


T  -  (2  )  T  t*(T-t* )  >  T  {  1  -  (_2)  T2  }  ^  0 

t2 


(4.22) 


recalling  that  0  t*  '  T.  Hence, 


d  t*(T)  -  0  for  T  r  [0,  Ts  ] 
"dT  v"r 


(4.23) 


A  - 

*  v  •/ 
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FIGURE  3:  COMPARISON  BETWEEN  Tr,  T*  AND  T* 


Now  if  i  1,  since  d  t*(T)  i-  0,  we  can  in  fact  deduce  that 

dT 


d  t*(T)  0  for  T  [0,  T r  (,»  )],  a  -  1 


and  hence  the  value  of  T  which  maximizes  t*(T)  is  Tr  (i  ). 

We  will  make  one  more  observation  before  presenting  a  complete  solution 
Recall  that  Tr(  )  for  the  case  <  ■  1  can  be  found  by  solving  (4.13)  and 
(4. 15) ,  that  is, 


(1-  )  -  1  (Ts/Tr)2 


exp  — 


1  (.  1  -  ) 


exp  ) 

(1-  ) 


Now,  let  (■;  =  1  and  s  =  1- 
equations 


(l-s)s  =  1  (Ts/Tr) 2 


in  (4.25)  and  (4.26),  to  obtain  the  pair  of 


(-■-!) 

s(l-s) 


81  (s) 


But  1,  hence  (4.27)  and  (4,28)  are  solved  by  Tr  (  t)  *  Tr  ( ^  ) 

=  Tj.il/.)  from  our  :,rfv  ions  work. 


Now,  in  order  to  find  the  bias  t*  induced  by  the  delay  Tr(a  ),  we  need 
to  solve  equation  (4.14)  with  T  =  Tr(a  ).  We  note  that  there  must  be  two 
values  of  which  solve  (4.14),  i.e.  !■'  =  g32  (a  ),  which  corresponds  to 
the  "flat  spot”  which  defines  Tr(  a),  and  the  value  of  F  which  corresponds 
to  the  maximum  bias  t*.  Define  the  function  g3  by 


>  =  1  -  exp  2(Tr/Ts)2 


)  =  g3(Tr;^) 


and  define  g4(  a)  by 


g4^  t) 


=  g3(Tr(  O;'-' 


f  gi  (•"’ 


)} 


( 


Then  g4  is  defined  and  single-valued  for  all  at  <  1.  Figure  4  shows  a  graph 
of  g3  for  representative  values  of  (Tr/Ts)  superimposed  on  a  graph  of 
g]_.  Figure  5  shows  a  graph  of  g4* 


Hence  from  (4.12)  we  conclude 


t*  =  g4(  O  Trd  )  *  g 4(  ‘  )  Tj-d/ot) 


( 


We  have  now  obtained  a  complete  solution  for  the  problem  of  choosing  a 

delay  so  as  to  maximize  the  bias  induced  in  a  victim  radar.  The  solution 
divides  naturally  into  two  cases: 


(i)  ■  -  Es/Ed  1 


(ii)  :  =  Eg/ Ed  1 


Optimal  bias  =  topt  =  Ts  g31(  0 
Optimal  delay  -  TQp£  -  t0p£  +  Ts/»2 

=  Ts  (g^d)  +  i  )  ( 

v  2 

Optimal  delay  =  Topt  =  1  Ts  [  1  ]1/2 

v'2  ~~ 


where 

Optimal  bias 


g].1  (  t)  or  ;  =  gjid/u) 
topt  ~  S4d)  ^opt 


Graphs  of  the  solutions  (T0  pt^^g)  an<^  d0pt/Ts)  are  shown  in  Figure  6. 


^  S 


LU 


94  (Es/Ed) 
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FIGURE  5:  GRAPH  OF  g3  FOR  SEVERAL  VALUES  OF  Tr/T 
SUPERIMPOSED  ON  GRAPH  OF  g[ 


i'r  v--  *-~  * 1  r'-r '  '*■■■*•* 


6.0  CONCLUSIONS 


We  set  out  to  derive  some  basic  "rules  of  thumb"  regarding  the  possible 
performance  and  response  time  of  a  hypothetical  ECM  system  which  used  the 
range-gate  (or  velocity-gate)  pulloff  countermeasure.  Our  analysis  has  been 
based  on  a  standard  statistical  model  for  the  radar  return  from  a  slowly 
fluctuating  point  target  in  the  presence  of  white  noise.  Using  our  model, 
which  also  included  a  deception  signal  representing  the  range-gate  pulloff 
ECM,  we  derived  an  expression  for  the  expected  power  envelope  of  the  filtered 
combined  signal.  We  saw  that  target  locations  could  be  identified  as  local 
maximum  of  this  power  envelope.  Moreover,  we  showed  that  the  same  basic  model 
was  also  appropriate  for  the  analysis  of  velocity-gate  pulloff. 

We  then  explored  our  model  in  more  detail  for  the  specific  case  in  which 
the  transmitted  signal  could  be  represented  by  a  gaussian  pulse.  We  saw  that 
small  deception  delays  resulted  in  biased  estimates  of  target  position; 
however,  as  the  delay  was  increased,  a  threshold  was  reached  beyond  which  the 
victim  radar  could  resolve  the  jammer  as  a  separate  target.  Now,  in  order  for 
the  countermeasure  to  be  successful,  one  must  be  able  to  move  the  apparent 
target  position  sufficiently  far  from  the  true  position  that  the  victim  radar 
cannot  keep  both  the  true  and  biased  positions  within  its  range-gate.  As  one 
might  have  expected,  the  determining  factors  as  to  whether  this  can  be 
accomplished  appear  to  be  the  pulse  width  of  the  transmitted  signal,  and  the 
signal-to-jammer  power  ratio  (calculated  at  the  receiver). 

A  number  of  possibilities  exist  for  extending  or  enhancing  this  model. 

For  example,  at  the  moment  the  additive  deception  signal  is  restricted  to  an 
amplified  and  delayed  version  of  the  signal  transmitted  by  the  radar.  One 
could  allow  a  more  general  type  of  deception,  such  as  one  which  included 
variations  in  the  pulse  width.  Indeed,  one  could  go  so  far  as  to  consider  an 
adaptive  jammer  which  constructs  a  jamming  signal  (based  on  the  measured 
transmitted  waveform)  designed  to  induce  the  maximum  deception  in  the  victim 
radar.  Needless  to  say,  we  suspect  this  problem  would  prove  to  be  very 
difficult.  Another  direction  might  be  to  attempt  a  simultaneous  solution  of 
the  range-gate  and  velocity-gate  problems;  reformulating  the  model  in  terms  of 
the  ambiguity  function  would  be  a  possible  approach.  One  might  also  wish  to 
add  a  more  realistic  description  of  the  dynamics  of  the  victim  radar,  or 
include  a  model  of  the  observer  needed  by  the  ECM  system  to  measure  the 
parameters  of  the  transmitted  signal. 

We  caution  the  reader  once  again  that  this  paper  is  aimed  primarily  at 
those  who  want  to  do  a  "top  level"  analysis  of  the  effectiveness  of  the 
range-gate  pulloff  ECM,  particularly  simulation  designers.  We  should  also 
point  out  that  there  are  some  obvious  limitations  which  must  be  considered 
before  applying  these  results  to  specific  problems.  For  example,  this  model 
would  certainly  not  be  an  appropriate  one  for  describing  a  tracking  gate  which 
used  sophisticated  signal  processing.  Also,  it  is  important  to  bear  in  mind 
that  this  model  describes  signals  by  their  statistical  averages;  consequently, 
one  should  not  attempt  to  draw  conclusions  about  the  pulse-by-pulse  behaviour 
of  such  systems  from  this  analysis.  However,  used  carefully  in  the  context 
for  which  It  was  intended,  we  believe  that  our  model  can  provide  a  useful  tool 
for  the  F.W  systems  analyst. 
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APPENDIX  A 
COMPUTER  I 'OGRAMS 


The  numerical  results  used  in  this  paper  (which  appear  primarily  in  the 
graphs)  were  generated  by  programs  written  in  the  APL  programming  language 
using  a  DEC-2020  minicomputer.  In  this  appendix  we  will  briefly  describe  each 
of  these  programs;  these  descriptions  will  be  followed  by  a  comprehensive 
listing . 

(a)  CHANGE:  CHANGE  accepts  a  boolean  vector  (i.e.  one  consisting  of 
0's  and  1' s) ,  which  we  will  denote  by  (a^,  a2,....an),  as 

an  argument.  It  returns  a  vector  of  indices  i  corresponding  to 
those  a^  which  satisfy  a^  t  a-j^. 

(b)  INV :  INV  is  an  infix  function  which  takes  as  its  first  argument 
the  name  of  the  function  to  be  inverted.  The  second  argument  is 
a  three  element  array  composed  of  the  value  at  which  the  inverse 
function  is  to  be  evaluated  followed  by  the  interval  over  which 
the  inverse  is  to  be  calculated.  That  is, 

f  INV  (x,a,b)  =  f_1(x)  for  x  f([a,b]) 

=  </,  for  x  i  f ( [ a , b] ) 

Note  that  f”-*-(x)  may  be  set-valued.  The  basic  algorithm  used 
to  compute  the  inverse  function  is  quite  simple.  The  interval 
[a,b]  is  subdivided  into  a  number  of  intervals  [a^,  b^ ] ,  and 
the  program  checks  to  see  if  x  f([aj,  b^]).  If  It  Is,  the 
process  is  repeated  until  the  size  of  the  subinterva 1( s)  is  less 
than  a  pre-selected  tolerance. 

(c)  G1  :  G1  is  the  function  defined  by  (4.15),  i.e. 

81  =  1~  '  exP  (_  -  .0 

Y-'*~ 

Since  g^  has  asymptotes  at  0  and  1,  traps  are  set  to  catch 
these  values. 

(d)  G2:  G2  is  the  function  defined  by  (4.20),  i.e. 

g2  ( x)  =1  exp  • x^  -  1 

/  ,  1 

1  w  X  — 

Since  g2  has  an  asymptote  at  0,  a  trap  is  set  to  catch  small 
values  of  x. 


G3:  G3  is  an  infix  representation  of  the  function  defined  by 
(4.29),  i.e. 

(Tr/Ts)  G3  F  =  g3(Tr ; '  ) 


g3(Tr;  F)  =  1  -F  exp  { 2(Tr/Ts)2  (F  -  1)} 

~y~  2 

Note  that  the  asymptote  at  0  is  trapped. 

G4:  G4  is  the  function  defined  by  (4.30),  i.e. 


g4  (J)  =  {  .’  |  1  =  g3(Tr(  >);/’)  ,  a  £  gl  (")} 

G4  can  accept  vector-valued  arguments;  however,  since  G4  works  by 
invoking  INV,  the  input  vector  is  disassembled,  the  inverse  is 
calculated  for  each  scalar  that  results,  and  the  answer  is 
assembled  into  an  output  vector.  Note  that  any  function  which 
invokes  INV  is  accompanied  by  a  SUSPEND  block  which  returns 
control  to  the  user  if  INV  returns  the  empty  set  or  multiple 
values.  Often  this  happens  because  of  numerical  considerations 
even  when  the  inverse  function  exists  and  is  single-valued. 

Rather  than  provide  complicated  code  to  deal  with  such  problems, 
it  is  easier  when  using  an  interpreter  like  APL  to  return  control 
to  the  user  and  let  him  sort  out  the  difficulty. 

SOLUTION:  SOLUTION  finds  the  values  of  T0pt  and  t0pt  (as 
defined  by  (4.32))  which  correspond  to  various  values  of 
(Eg/Efj)  ranging  from  -50  to  50  dB.  LOOPA  solves  part  (i)  of 
(4.32),  and  LOOPB  solves  part  (ii).  Note  the  presence  of  the 
SUSPEND  block  once  again  since  INV  Is  invoked  frequently. 


V  I NDEX*CHANGE  MASK  ;I  ;J  ;S AVEMASK  ;F LAG 

[  1  ]  A 

[2]  a  MASK  IS  A  VECTOR  OF  O'S  AND  l'S 

[3]  a  CHANGE  FINDS  INDICES  CONTAINING  VALUES  WHICH 

[4]  A  DIFFER  FROM  PRECEDING  VALUE 

[5]  A 

[  6  ]  INDEX*  i  0  ;  MAS  K+-MAS  K  ,  0  -,J*Q  -,FLAG*Q  ;S  AVEMASK*MASK 

[  7 ]  LOOP :*(0 -qMASK*( I *MASK 1 1 ) +MASK) /EXIT 

[8]  INDEX*INDEX,I+J+~1+INDEX 

[9]  *{0  =  pMASK*U*MASK\0)  +MASK)  /EXIT 

[10]  +LOOP 

[11]  EXIT :+(0*F LAG) /FINISH 

[12]  -J*  0  ;INDEX*\  Q  ;S  AVEIN*I  NDEX  ;MASK*(  ~S AVEMASK  )  ,  0  ; 

[13]  *LOOP 

[14]  FINISH :INDEX*(Q*INDEX) /INDEX+S AVEIN , INDEX 

[15]  INDEX *( ( oSA V  EM ASK ) xINDEX ) / INDEX*INDEX [ \INDEX ] 

[  16  ]  G 
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7  1+F’JNC  IN  V  V  ALU  E  ;LEFT  ;R IGHT  ;U N  ;X  ;Y  \MASK\I  NDEX  ;0K  ;  NOTOK 
[  1  ]  p 

[2]  p  IN T  FINDS  INVERSE  FUNCTION  EOF.  FUNC  EVA  LU  A  TEL 

[3]  fl  AT  FIRST  ELEMENT  OF  ARRAY  VALUE.  SECONL  AND 
[ 4  J  p  THIEL  ELEMENTS  OF  VALUE  DEFINE  INTERVAL 

[5]  p  OF  INTEREST. 

L  6  ]  p 

[7]  p  CHECH  TO  SEE  IF  VALUE  IS  LEGAL  FORMAT 

[8]  p 
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